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Abstract 

The preliminary data of the PHENIX cohaboration for the scaled variances of charged hadron 
multiplicity fluctuations in Au+Au at ^/s = 200 GeV are analyzed within the model of independent 
sources. We use the HSD transport model to calculate the participant number fluctuations and 
the number of charged hadrons per nucleon participant in different centrality bins. This combined 
picture leads to a good agreement with the PHENIX data and suggests that the measured multi- 
plicity fluctuations result dominantly from participant number fluctuations. The role of centrality 
selection and acceptance is discussed separately. 
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The event-by-event fluctuations in high energy nucleus-nucleus (A+A) collisions (see e.g., 



the reviews 
(see e.g., Refs. 



, 1^ are expected to provide signals of the transition between different phases 
I) and the QCD critical point [5|]. In the present letter we study the 
charged multiplicity fluctuations in Au+Au collisions at RHIC energies. The preliminary 
data of the PHENIX collaboration [6] at y/s = 200 GeV are analyzed within the model 
of independent sources while employing the microscopic Hadron- String-Dynamics (HSD) 
transport model [3, B] to define the centrality selection and to calculate the properties of 
hadron production sources. 

The centrality selection is an important aspect of fluctuation studies in A+A collisions. At 
the SPS fixed target experiments the samples of collisions with a fixed number of projectile 
participants N^"-' can be selected to minimize the participant number fluctuations in the 
sample of collision events. This selection is possible due to a measurement of the number of 
nucleon spectators from the projectile, Ng"°-' , in each individual collision by a calorimeter 
which covers the projectile fragmentation domain. However, even in the sample with Np "'' = 
const the number of target participants fluctuates considerably. In the following the variance, 
Var{n) = (n^) — (n)^, and scaled variance, u = Var{n)/ (n) , where n stands for a given 
random variable and (■ ■ ■ ) for event-by-event averaging, will be used to quantify fluctuations. 
In each sample with iVp = const the number of target participants fluctuates around its 
mean value, {N^p"^^) = N^°-^ , with the scaled variance cup^^. Within the HSD and UrQMD 
transport models it was found in Ref. Isi] that the fluctuations of N^p^^ strongly influence 
the charged hadron fluctuations. The constant values of N^°-^ and fluctuations of N*p^^ lead 
also to an asymmetry between the fluctuations in the projectile and target hemispheres. The 
consequences of this asymmetry depend on the A+A dynamics as discussed in Ref. 10 1. 

In Au+Au collisions at RHIC a different centrality selection is used. There are two kinds 
of detectors which define the centrality of Au+Au collision, Beam-Beam Counters (BBC) 
and Zero Degree Calorimeters (ZDC). At the cm. pair energy ^/s = 200 GeV, the BBC 
measure the charged particle multiplicity in the pseudorapidity range 3.0 < |?7| < 3.9, and 
the ZDC - the number of neutrons with |?7| > 6.0 ^. These neutrons are part of the nucleon 
spectators. Due to technical reasons the neutron spectators can be only detected by the ZDC 
(not protons and nuclear fragments), but in both hemispheres. The BBC distribution will be 
used in the HSD calculations to divide Au+Au collision events into 5% centrality samples. 
HSD does not specify different spectator groups - neutrons, protons, and nuclear fragments 
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such that we can not use the ZDC information. In Fig. [T] (left) the HSD results for the BBC 
distribution and centrality classes in Au+Au collisions at ^/s =200 GeV are shown. We 
find a good agreement between the HSD shape of the BBC distribution and the PHENIX 
data [6]. The experimental estimates of (Np) for different centrality classes are based on 
the Glauber model. These estimates vary by less than 0.5% depending on the shape of the 
cut in the ZDC/BBC space or whether the BBC alone is used as a centrality measure ^. 
Note, however, that the HSD (Np) numbers are not exactly equal to the PHENIX values. 
It is also not obvious that different definitions for the 5% centrality classes give the same 
values of the scaled variance ujp for the participant number fluctuations. 
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FIG. 1: HSD model results for Au+Au collisions at y/s = 200 GeV. Left: Centrality classes defined 
via the BBC distribution. Right: The average number of participants, {Np), and the scaled variance 
of the participant number fluctuations, ujp, calculated for the 5% BBC centrality classes. 



Defining the centrality selection via the HSD transport model (which is similar to the 
BBC in the PHENIX experiment) we calculate the mean number of nucleon participants, 
(Np), and the scaled variance of its fluctuations, up, in each 5% centrality sample. The 
results are shown in Fig. [H right. The Fig. [2] (left) shows the HSD results for the mean 
number of charged hadrons per nucleon participant, rii = {Ni)/{Np), where the index i 
stands for "— ", "+", and "ch", i.e negatively, positively, and all charged final hadrons. 
Note that the centrality dependence of Ui is opposite to that of up: rii increases with {Np), 
whereas up decreases. 

The PHENIX detector accepts charged particles in a small region of the phase space with 
pseudorapidity |?7| < 0.26 and azimuthal angle (p < 245° and the pt range from 0.2 to 2.0 
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GeV/c j^. The fraction of the accepted particles qi = {Nf^^) / {Ni) calculated within the 
HSD model is shown in the r.h.s. of Fig. [2l According to the HSD results only 3 ~- 3.5% of 
charged particles are accepted by the mid-rapidity PHENIX detector. 



15 



10 




0.035 



0.030 



100 



200 
<Np> 



300 



400 




FIG. 2: HSD results for different BBC centrality classes in Au+Au collisions at ^/s = 200 GeV. 
Left: The mean number of charged hadrons per participant, = (Ni) / {Np) . Right: The fraction 
of accepted particles, qi = {N°-^^) / {Ni) . 

To estimate the role of the participant number event-by-event fluctuations we use the 
model of independent sources (see e.g., Refs 



ifli) 



UJi 



UJ, + Hi LUp . 



(1) 



The first term in the r.h.s. of Eq. ([T]) corresponds to the fluctuations of the hadron mul- 
tiplicity from one source, and the second term, up, gives additional fluctuations due to 
the fluctuations of the number of sources. As usually, we have assumed that the number of 
sources is proportional to the number of nucleon participants. The value of Ui in Eq. ([T]) is 
then the average number of i'th particles per participant, rii = {Ni)/{Np). We also assume 
that nucleon-nucleon collisions define the fluctuations u* from a single source. To calculate 
the fluctuations cuf^ '^ in the PHENIX acceptance we use the acceptance scaling formula (see 
e.g., Refs. [l|,y,Q): 



a;""" = I - qi + qiUJi 



(2) 



where qi is the fraction of the accepted i'th hadrons by the PHENIX detector. Using Eq. ([T]) 
for LJi one finds, 



- Qi +qi ^*i + Qi rii Up . 



(3) 
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The HSD results for up (Fig. [H right), rii (Fig. [21 left), qi (Fig. [2l right), together with the 
HSD nucleon-nucleon values, u'L = 3.0, uj\_ = 2.7, and tu*^ = 5.7 at i/i = 200 GeV, define 
completely the results for uj"''^'^ according to Eq. ([3]). We find a surprisingly good agreement 
of the results given by Eq. (j3]) with the PHENIX data shown in Fig. [3l Note that the 
centrality dependence of tuf'^^ stems from the product, rii ■ ujp, in the last term of the r.h.s. 
of Eq. ([3D. 
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FIG. 3: The scaled variance of charged particle fluctuations in Au+Au collisions at ^/s = 200 GeV 
with the PHENIX acceptance. The circles are the PHENIX data while the open points (con- 
nected by the solid line) correspond to Eq. ([3]) with the HSD results for ujp, rii, and qi. 

In summary, the preliminary PHENIX data for the scaled variances of charged hadron 
multiplicity fluctuations in Au+Au collisions at -^i = 200 GeV have been analyzed within 
the model of independent sources. Assuming that the number of hadron sources are propor- 
tional to the number of nucleon participants, the HSD transport model was used to calculate 
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the scaled variance of participant number fluctuations, uip, and the number of i'th hadrons 
per nucleon accepted by the mid-rapidity PHENIX detector, qiUi, in different Beam-Beam 
Counter centrality classes. The HSD model for nucleon-nucleon collisions was also used to 
estimate the fluctuations from a single source, uj*. We find that this model description is in 
a good agreement with the PHENIX data [6]. In different (5%) centrality classes up goes 
down and g^nj goes up with increasing (Np). This results in non-monotonic dependence of 
^acc ^jY^^ gggj^ ^Yie PHENIX data. 

We conclude that both qualitative and quantitative features of the centrality dependence 
of the fluctuations seen in the present PHENIX data are the consequences of participant 
number fluctuations. To avoid a dominance of the participant number fluctuations one needs 
to analyze most central collisions with a much more rigid (< 1%) centrality selection. The 
statistical model then predicts uj± < 1 [ill], whereas the HSD transport model predicts the 
values of uj± much larger than unity at ^/s = 200 GeV [l2|. To allow for a clear distinction 
between these predictions it is mandatory to enlarge the acceptance of the mid-rapidity 
detector up to about 10% (see the discussion in Ref. |l2|). 
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